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Abstract 

We consider the semiclassical theory in a joint phase space of spin and orbital 
degrees of freedom. The method is developed from the path integral using the spin- 
coherent-state representation, and yields the trace formula for the density of states. 
We discuss the limits of weak and strong spin-orbit coupling and relate the present 
theory to the earlier approaches. 



1 Introduction 

The subject of spin-orbit interaction in solid state physics has been attracting much 
attention recently due to potential applications in spin-based electronic devices 
IT2] . It has often proved advantageous, on the other hand, to exploit the semiclassical 
description of electrons in mesoscopic systems in the ballistic regime |28| . Hence there 
is a substantial interest in extending the semiclassical theory to include spin. Several 
attempts were made in the last decade in this direction. Littlejohn and Flynn 
revised and improved the asymptotic theory of coupled wave equations and applied 
it to systems with standard spin-orbit coupling. Their approach, however, relies on a 
finite strength of the spin-orbit interaction, and thereby becomes invalid at the points 
in classical phase space where the interaction vanishes — the so-called mode-conversion 
points. Frisk and Guhr ^3] found that in certain cases this problem can be corrected 
by a heuristic procedure. Bolte and Keppeler [2] studied the opposite situation when 
the interaction is weak. They derived the trace formula for the Dirac equation and 
Pauli equation with spin 1/2. Based on the assumption that the orbital motion is not 
affected by spin, their approach was used to explain the anomalous magneto-oscillations 
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in quasi- two-dimensional systems (201 • A comparative analysis of the above-mentioned 
semiclassical methods and their applications to specific systems can be found in |2j. 

In a new semiclassical approach to spin-orbit coupling }27| the orbital and spin 
degrees of freedom are combined in a joint (extended) phase space. Then the whole 
arsenal of the semiclassical methods developed for spinless systems can be applied 
there. The description of spin by continuous variables is achieved by using the basis 
of coherent states. After a path integral containing the spin and orbital coordinates is 
constructed, it can be evaluated in the stationary-phase approximation, yielding the 
semiclassical propagator or its trace. 

A detailed investigation of the method of 27 a , with the focus on the trace formula 
for the density of states, is the subject of the present paper. We will argue that the 
earlier approaches can be put under one roof by this new theory. (A step in this 
direction was made in [HE]-) Our discussion will be limited to Hamiltonians linear in 
spin. For such systems, as we will show, the semiclassical treatment is valid not only 
in the large-spin limit, but also for a finite spin, including spin 1/2. 

We expect that the semiclassical method presented below can be effectively used in 
systems with relatively strong spin-orbit interaction, such as p-InAs or InGaAs-InAlAs 
heterostructures or atomic nuclei [I]. Another field of application is molecular 
dynamics, where similar methods were developed to map the discrete electronic states 
to continuous variables and then to treat those semiclassically |33j . 

Our paper is organized as follows. In the introductory section [21 we briefly review 
the derivation of the trace formula without spin, starting from the path-integral rep- 
resentation for the propagator. Recently proposed by Sugita [52] . this new procedure 
reproducing the well-known result of Gutzwiller j£j will later be generalized to include 
spin. In Sec. |3] we define the spin coherent states and write the quantum-mechanical 
partition function in the path-integral representation for a system with spin-orbit cou- 
pling. Sec. El is the main part of this work, where three versions of the trace formula 
are derived. In the first two instances we integrate over the spin variables in the path 
integral exactly and then apply the stationary-phase approximation. Consequently, for 
weak spin-orbit interaction fSec. I4.l1 and Appendix ^J) we obtain the generalization of 
the Bolte-Keppeler trace formula to arbitrary spin, and for strong coupling (Sec. I4.2JI 
we recover, on a restricted basis, the result by Littlejohn and Flynn (at the moment we 
cannot recover their "no-name" phase). In the third case (Sec. 14.31) the path integral 
is evaluated by the stationary-phase method in both spin and orbital variables. Then, 
using Sugita's approach, we express the oscillating part of the density of states as a sum 
over the periodic orbits in the extended phase space. Special attention is paid to the 
Solari-Kochetov phase correction |3U1 1221 134j arising in the semiclassical limit of spin- 
coherent-state path integrals. In Sec. 14.41 we classify the regimes of weak, strong, and 
intermediate coupling, and show that the trace formula in the extended phase space 
works in all of them. (Appendix iBl illustrates explicitly the strong-coupling limit.) We 
also discuss and relax the large-spin applicability condition for the semiclassical ap- 
proximation. As an application of the general theory, in Sec. El we study two systems 
with Rashba-type spin-orbit interaction [S] and the Jaynes-Cummings model |19j . 
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2 Trace formula without spin 



In this section we summarize some basic facts related to the semiclassical trace formula 
and present its derivation for systems without spin We consider a system described 
by the Schrodinger equation 

Hi, k = E k ^ k (1) 

with a discrete energy spectrum {E k }- Its density of states g(E) = 8{E — E k ) can 
be subdivided into smooth and oscillating parts, i.e., 

g{E) = g(E) + 5g(E) . (2) 

From the semiclassical point of view, the smooth part g(E) is given by the contribution 
of all orbits with zero length and can be evaluated by the (extended) Thomas-Fermi 
theory [3]. Numerically it can be extracted by a Strutinsky averaging of the quantum 
spectrum jSJ. 

In this paper we assume the smooth part to be known and will be interested in the 
oscillating part 5g(E), which is semiclassically approximated by the trace formula 15 

5g sc (E) = ^Apo(E) cos ( ^S po (E) - ^ P oj • (3) 

po ^ ' 

The sum here is over all classical periodic orbits (po), including all repetitions of each 
primitive periodic orbit (ppo). S po (E) is the action integral and a po is the Maslov 
index of a periodic orbit. The amplitude A po (E) depends on the integrability and the 
continuous symmetries of the system. When all periodic orbits are isolated in phase 
space, the amplitude is given by JH] 

A P o(E) = -L Z VVO (4) 

Y / |det(M po -/ 2(d _ 1) )| 

where T ppo is the period of the primitive orbit, M po is the stability matrix of the periodic 
orbit, and d is the number of degrees of freedom. / is the unit matrix, the subscript 
denotes the dimensionality of the space where / acts. For two-dimensional systems 



det(M po — l2)\ = 2| sin(A po /2)|, where A po is called the stability angle of a periodic 
orbit. 

Recently Sugita [22] has proposed a re-derivation of Gutzwiller's trace formula di- 
rectly from the quantum partition function, avoiding the calculation of the semiclassical 
propagator. Since we utilize his approach in the following sections for systems with 
spin, we will briefly review it here. Sugita's starting point is the quantum partition 
function (or the trace of the quantum propagator) in the path integral representation 

Z(T) = j V[ ^}y ] ex P {^[q,p;T]/a}. (5) 

The path integral is calculated along closed paths (q(i),p(i)) in the 2d-dimensional 
phase space over a time interval T: 

mm s lim A A (M&M) , tl = a in. (6 ) 

(2nKY JV-oo-l-llll 2irh ) ' 11 K 1 

y ' 1=1 j=l V 7 
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In © 



ft[q,p;T] 



(p q-q p) - H(q,p) 



(7) 



is the Hamilton principal action function. W(q, p) is the classical Hamiltonian, i.e., the 
Wigner-Weyl symbol of H. The action form p • dq can be antisymmetrized along the 
closed paths. The Fourier-Laplace transform of Z(T) yields, after taking the imaginary 
part, the density of states 



9(E) 



1 



Im 



7T 



1 



e iET/h Z (T)dT 



(8) 



The path integral (J5J) receives its largest contributions from the neighborhoods of 
the classical paths, along which the principal function 1Z is stationary according to 
Hamilton's variational principle 51Z = 0. The first variation hereby yields the classical 
equations of motion 

V--— O) 

One may evaluate the integrations in © using the stationary-phase approximation, 
which becomes asymptotically exact in the classical limit 1Z 3> H. The semiclassical 
approximation of the partition function Z(T) then turns into a sum over all classical 
periodic orbits with fixed period T 



ZsciT) 



^2 e m P° /h 

po 



^ exp {in p 2 J[r,;T]} 



(10) 



where lZ po are the principal functions (J7J), evaluated now along the classical orbits. The 
functional of the second variations is 



K p 2 J[r,,T] = £ 



dt 



(A,p) 



(5q, 5p) I \/~h and J 



is the 2d-dimensional unit symplectic matrix. H,^ 



where rj is the 2d-dimensional phase-space vector of small variations r) 

h 
-h 

is the second variation of the classical Hamiltonian Ti, calculated along the periodic 
orbits: 

d 2 n 

piPj- 



1 d 



^ j dqidq 



d 2 n nx d 2 n 



dqidpj 



dpidpj 



(12) 



Note that Z SC (T) does not include the contribution of the zero-length orbits. 

After a stationary-phase evaluation of the Fourier-Laplace integral (jSJ with Z SC (T) 
instead of Z(T), one finally obtains Gutzwiller's trace formula © where the actions 
Sp (E) = lZp +ET po are calculated at fixed energy E and the periods of the orbits are 
dSp /dE. The monodromy matrix M po is defined by rj(T po ) = M po in terms 



L po 



of the solutions of the linearized equations of motion r] = J dTi^'/di], which are purely 
classical. After removing the trivial parabolic block from M po that appears due to the 
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time translation symmetry, one obtains the reduced 2{d — l)-dimensional monodromy 
matrix M po . The latter enters formula (j3J) and since it contains information about the 
stability of periodic orbits [§] it is often referred to as the stability matrix. For the 
Maslov indices a po Sugita has also given general formulae |32| . 

3 Spin coherent states and quantum partition 
function 

We intend to generalize the path-integral representation for the partition function (J5J) 
to include the spin degree of freedom. The main issue then is to be able to describe 
spin on the quantum-mechanical level by a continuous variable. One can achieve this 
by using the overcomplete basis of spin coherent states. The path integral for a system 
with spin in the SU(2) spin-coherent-state representation originally appeared in a paper 
by Klauder as an integral on the sphere §> 2 . Kuratsuji et al. 23, have represented 
it as an integral over paths in the extended complex plane C . 
The SU(2) coherent state \z;s) for spin s is defined by 

\z; s) = (1 + |z| 2 ) _s exp(zs + ) \s, m s = — s) , (13) 
s-\s, m s = — s) = , (14) 

where z G C 1 is a complex number. The spin operators s± = §i ± i§2 and S3 are the 
generators of the spin su(2) algebra: 

[s3,s±]=±s ± , [s + ,s_]=2s 3 . (15) 

\s,m s ) are the eigenstates of S3. From the group-theoretical point of view, s S N/2 
labels irreducible representations of SU(2). 

The irreducibility, as well as the existence of the SU(2) invariant measure d/i s , 
ensures that the resolution of unity holds in the spin-coherent-state basis: 

J \z; s)(z;s\ dfi s (z*,z) = J 2s+ i , 
dfj, s (z*,z) = ,^ S + , \, 2 d 2 z . (16) 

7T (1 + \Z\ Z ) Z 

This turns out to be the most important property of spin coherent states that allows 
for the path-integral construction. The measure d(i s takes account of the curvature of 
the sphere § 2 . In what follows, we denote \z;s) simply by \z). 

Let us now consider a quantum Hamiltonian with spin-orbit interaction 

H = H + H s = H (<i,p) + hs-C(ei,p), (17) 

where s = (si,S2, S3) and C = (Ci, 62,63) is a vector function of the coordinate and 
momentum operators q, p. Thereby we assume the most general form of spin-orbit 
interaction linear in spin. We write the expression for the respective quantum propaga- 
tor in terms of a path integral in both the orbital variables q, p and the spin-coherent- 
state variables z,z*. Imposing periodic boundary conditions on the propagation and 
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thus integrating over closed paths, we arrive at the expression for the partition func- 
tion [cf. ©]: 

Z(T) = J V[ ^y V» a [z\z] exp{iR,[<i,p,z*,z;T}/h}. (18) 
The Hamilton principal action function now includes the symplectic 1-form due to spin: 



K[q,p,z*,z;T] =j 



1 zz* — z* z 

-(p-q-q-p) + i(i + \ z \2j ~ H (^Pi z *^ z ) 



dt , (19) 



and the path integration in (|18|) is taken over the 2{d+ l)-dimensional extended phase 
space: 

P[q]P[p] -r, r* 1 v 11 T I dqj(U)dpj(ti) j 

V ' 1=1 3=1 

where the time interval T is divided into N time steps t\ = IT/N . 

The c-number Hamiltonian W(q, p, z* , z) appearing in the integrand of (|19|) is 

H(q, p, z*, z) = H (q, p) + W s (q, p, z*, z) = Ho(q, p) + fen(z*, z) • C(q, p) , (21) 

where TCo(q, p) and C(q, p) are the Wigner-Weyl symbols of the operators ifo and C, 
and n = (ni, 712,723) = (z\s\z)/s is the unit vector of dimensionless classical spin. Its 
components in terms of z* and z are given by 

ni + in 2 = 2z*/(l + |z| 2 ), n 3 = -(1 - |^| 2 )/(1 + k| 2 ). (22) 

One may recognize the stereographic projection from a point n on the unit Bloch sphere 
onto a point (Re z, — Im z) on the complex plane, whereby the south pole is mapped 
onto its origin z = 0. 

We remark that the choice of the c-number Hamiltonian is not unique and, in 
the discretized path integral, is linked to a specific discretization prescription |24| . 
For symbolic manipulations in the continuous limit, it is often convenient to use the 
Wigner-Weyl symbol of the quantum Hamiltonian as the c-number Hamiltonian. Note, 
however, that sn is the covariant symbol of the spin operator s, and not the Wigner- 
Weyl symbol. The consequences of the present choice of H(q, p, z*,z) (f2*T|) for the 
semiclassical theory are discussed in Sec. 14.41 



4 Semiclassics with spin-orbit coupling 

In this section we consider the semiclassical approach for systems with spin-orbit in- 
teraction. We begin with the limit of weak coupling (Sec. 14. 1|) . calculating the spin 
part of the path integral exactly. The trace formula we obtain is a generalization of 
the results by Bolte and Keppeler [3] to arbitrary spin. A similar approach leads, un- 
der certain restrictions, to the trace formula in the adiabatic (strong-coupling) limit 
(Sec. 14.2(1 . that was studied by Littlejohn and Flynn |25j . In Sec. 14.31 we present a 
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new method that treats both orbital and spin degrees of freedom semiclassically |27j . 
With the latter approach one can go beyond the limits of weak and strong coupling. 
The relationship between the methods, in connection with the coupling strength, is 
analyzed in Sec. 14.41 

4.1 Weak-coupling limit: quantum- mechanical treatment 
of spin 

We derive here the semiclassical trace formula for the Hamiltonian (|17jl in the asymp- 
totic limit h — > 0. In addition to the standard semiclassical requirement IZ/h S> 1, 
which is automatically fulfilled, this limit also implies that the spin-orbit interaction 
energy is small: 

|w.(q,P,*V)|<|w (q,p)l, ( 23 ) 

i.e., the system is in the weak-coupling regime. 

It is convenient to separate the terms that explicitly contain h in the Hamilton 
principal action function (|19|h Representing 

TZ[q,p,z*,z;T}=n [c L ,p;T}+h S TZ 1 [c L ,p,z*,z;T}, (24) 

where the unperturbed part TZq is given by ((7J) with the Hamiltonian 7io(c[,p) an d 

rT 



TZ 1 [c l ,p,z*,z;T\ = j> 



n(z*,z) ■ C(q,p) 



dt , (25) 



with 



i{l + \z\ 2 ' 

we can write the partition function (j!8j) in the form 

Z(T) = J M,[cl, p; T] exp {i^ [q, p; T]/h} (26) 

M s [c L ,p;T} = Jv^iz^zjexpiisTZ^p^^z^T}}. (27) 

Now we can evaluate Z(T) in the stationary-phase approximation in q and p. Since 
s1Z\ <C TZo/h in the case under consideration, -M s [q, p; T] is a slowly varying functional 
of q and p and need not be taken into account when evaluating the stationary-phase 
condition. Applying the results of Sec. [2 we obtain the semiclassical approximation 
to the partition function as a sum over the unperturbed periodic orbits (determined 
entirely by Ho) 

Z SC (T) = £ M po {T) 4°) (T) exp {iK$/H - i$*/2) , (28) 



po 



where A p ° \ Vp°J , TZ^J = IZo [<lpo , Ppo ', T] , and the modulation factor M. po (T) = 
M s \ Qpo , Ppo ] T] are evaluated for the unperturbed trajectory. Performing the Fourier- 
Laplace transform Q (also in stationary-phase approximation) we arrive at the trace 
formula in the weak-coupling limit (WCL) 

^ CL = ^M po {E)A$(E) cos [iS$/h - ^tt/2} , (29) 

po 
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which differs from the unperturbed trace formula for TCq only by the presence of the 
modulation factor M po (E) = M po {T^{E)) with T^(E) = dS$ /dE. 

The modulation factor is determined in Appendix ^] by a direct evaluation of the 
path integral (|27j) . For a Hamiltonian linear in spin, the problem is effectively re- 
duced to the calculation of the spin partition function in the time-dependent external 
magnetic field C(q(t),p(t)) = (2Re/(i), -2Im/(i), 2A(t)), determined by the path 
(q(t),p(t)). Then for a periodic orbit one finds 



sin [(s + 1/2) rj(E )] 
sin[r?(E)/2] 



M P o(E) = L ll/ P w^ (30) 



with 

- T (0) 

V (E)= / dt{2A^(t)-[(f^(t)yz + f^(t)z*]}. (31) 



J o 

z*(t) and z(t) are to be found from the first-order differential equation 

z = -i[2A<®(t)z + / (0) (t) - (fW(t))*z% z(0) = z(T), (32) 

which is simply the precession equation h = x n, n(0) = n(T) for the classical 
spin vector n [Eq. (|22j) ]. This equation is angle-preserving, i.e., it rotates the Bloch 
sphere without deforming it. Thus, the orientations of the Bloch sphere at t = 
and t = T are related by a rotation about some axis. The angle of rotation is rj(E) 
(Appendix EJ) . The choice of initial condition n(0) along this axis corresponds to a 
periodic solution (cf. jHH])- Thus, (|32j) has two periodic solutions, whose rj(E) have 
opposite signs. Nevertheless, Ai po (E) is well defined. Note that i](E) is equal to the 
part TZi [q^J , p p °J , z*, z; T^] of the Hamilton principal action function calculated along 
the periodic z*{t) and z{t). Although we are able to give a classical interpretation to the 
ingredients of Ai po (E), no stationary-phase approximation was used in its derivation. 

In the absence of spin-orbit interaction or external magnetic tied one finds A4 po (E) = 
2s + 1, i.e., the unperturbed trace formula is multiplied by the spin degeneracy factor. 

For s = 1/2 the modulation factor M po (E) was derived in 5 by a different method, 
which, as ours, treats the spin degrees of freedom on the quantum-mechanical level. 

Clearly, when the weak-coupling condition IJ23JI breaks down, the trace formula (|29|) 
is no longer valid. Indeed, in this case 1Z\ is large, and A4 s [q, p;T] will influence the 
stationary-phase condition for the integral (|26|) , 

Note that the representation of terms in the WCL trace formula as a product of the 
unperturbed part and the modulation factor [Eq. I|29j) ] remains valid even when the 
Hamiltonian is nonlinear in spin. However, the simple expression for the modulation 
factor Eq. ()3U|) is specific to the Hamiltonian linear in spin. 

4.2 Adiabatic limit: strong coupling 

With an exact integration over the spin degrees of freedom in the path integral we 
can also derive the trace formula in the adiabatic limit. This limit is defined by the 
requirement 

|C(q,p)|»T-i, (33) 
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where T or b is the period of the orbital motion. Since |C| is the frequency of precession 
of the classical spin about the instantaneous magnetic field C, Eq. ()33|) means that 
the spin motion is much faster than the orbital motion. The results of this section 
will be most useful in the strong-coupling limit (SCL) \H S \ ~ \Ho\ or |C| ~ \Tto\/h, 
where (|3*3*|) is automatically satisfied. Formally speaking, the SCL can be stated as the 
double- limit h — > 0, |C| — > oo with fr\C\ = const (cf. 0). 

We start with the representation (|2(i)) for the partition function and write the pref- 
actor in the form (Appendix 

Ms[q,p;T] = — — f -— = > exp{im s ri q,p;T}. 34 

sm{r ? [q,p;T]/2} 

Then Z(T) becomes a sum over polarizations 

z(t)= £ f V ^}y exp Wq ' p; T]/h + imsV [q> p; T]} ' (35) 

where the path integrals will be calculated by the stationary phase. The functional 
7] [q, p; T] is given by [cf. (jS^)] 

r, [q, p; T] = f T dt {2A - [f*z + fz*}} (36) 
J o 

with / = |C| sm9ce~ l ^ c and A = (|C|/2) cos 6c in terms of the polar angles 9c, 4>c of 
the vector C. The trajectory z{t) is one of the two periodic solutions of the equation 

ri = Cxn, |n| = 1, (37) 

which can be solved approximately in the adiabatic limit. We look for the solu- 
tion in the form n = n(°) + nW, where n(°) || C and nW _L C. Assuming ~ 

n(°V(|c|r orb ) < n(0) , 

we obtain 

m C m n(°) x n(°) , N 

n = ± TcT ' = T — icj — ' (38) 

Then (|36|) gives us the phase 

m s i,[q lPi T]=-m s |C| dt + ipf s) + 0(|C| T^)" 1 . (39) 
JO 

Here the second term is the Berry phase 

fT 



r 1 

ipf s) =m s (1 + cos C ) <Pc dt. (40) 
Jo 

When the path integrals (|35[) for Z(T) are evaluated by stationary phase, the Berry 
phase if ^ ls ^ ~ 1 does not play a role in the determination of the stationary-phase point. 
In the SCL, i.e., when |C| ~ \7io\/h, the first term of (|3*9")) must be varied together 
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with TZo/h in order to derive the stationary-phase condition. Then after the standard 
procedure (Sec. |2j) we obtain the trace formula 

^ S S C CL = £ 6g£>\ (41) 

m s =— s 

where each polarized density of states 5gi™ 3 ^ is given by © with the classical dynamics 
controlled by the effective Hamiltonian 

s) (q, p) = H (q, P) + m s ^|C(q, p)| (42) 

and the Berry phase added. An independent semiclassical derivation of this result is 
given in Appendix iBl 

Actually, this result is valid only if C(q, p) depends on either qj or pj, but not 
both, for each j = l,...,d. In the latter case such symbolic manipulations with 
path integrals do not seem to be valid, and, as a result, the "no-name" term found 
by Littlejohn and Flynn 25 is missing. The problem of recovering this term in the 
path-integral approach has been known for some time and was already mentioned in 
Ref. |25j . Making some heuristic assumptions, Fukui ^1] reproduced the no-name term 
for the Jaynes-Cummings model in a framework similar to ours. It is not clear to us at 
this stage, whether his approach can be generalized to a wider class of Hamiltonians. 
We suppose that this term emerges due to the operator ordering, and hope to recover 
it making a careful analysis of the continuous limit of the path integral. 

It may happen that the adiabatic condition (|33j) is fulfilled in the WCL. The trace 
formula is well defined if the limit h\C\/TCo — > is taken before the limit |C| T or b — > oo. 
Hence, m s r] [q, p; T] does not contribute to the determination of the stationary point, 
and the orbital motion is governed entirely by TCq. The results of the previous section 
will be recovered with 17(E) given by ()39|) . i.e., the fast spin precession can be removed 
from r)(E) in the adiabatic limit. 

We note that expression (|39j) becomes invalid if the adiabatic condition Q33[) breaks 
down anywhere along the periodic orbit. In particular, this occurs in the mode- 
conversion points defined by C(q, p) = 0. 

4.3 Semiclassics in the extended phase space 

In the present approach, proposed recently in |27| . both the orbital and spin degrees of 
freedom are treated semiclassically. For the formal derivation we assume the standard 
semiclassical requirements of large action 1Z,S ^> h and large spin angular momen- 
tum, i.e., 

S = hs^h. (43) 

Later the second condition will be dropped. The spin-orbit coupling is allowed to be 
arbitrary now. Moreover, we show in Sec. 14.41 that the WCL trace formula (|29j) and 
the SCL trace formula (|4*Tj) can be reproduced within the approach of [2*7] . 

In Sec. [21 we presented the derivation of the spinless trace formula starting from the 
quantum partition function. Here we apply the same scheme to the spin-dependent 
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Z(T) of (|18() with the c-number Hamiltonian (|21j) responsible for the classical dynamics. 
At the end we should obtain the trace formula in its standard form (J2J). However, all 
its elements (periodic orbits, their actions, periods, and stabilities) are to be found 
from the dynamics in the extended phase space q, p, z*,z. In what follows we show 
how these ingredients can be determined. 

The first step is to obtain the classical equations of motion from the variational 
principle STZ[c[, p, z* , z] = 0, where, in particular, z* and z are varied independently. 
These generalized Hamilton equations are 



.(1 + 



|2\2 



&H 

dpi ' 

an 



Pi 



on 

dqi ' 

.(1 + 1* 



2\2 



on 



(44) 



(45) 



25 dz ' ' 25 dz* ' 

We emphasize that these equations describe the coupled dynamics of spin and orbital 
degrees of freedom and collectively determine the periodic orbits in the extended phase 
space. 

The second pair of equations are not redundant, as may appear. Actually, they 
prove the expected property that z* is indeed the complex conjugate of z along the 
classical paths. After we have established this, we can rewrite Eqs. (|45|) in terms of 
the real variables v = Im z* = — Im z and u = Re z* = Re z as 



. _ (l + \z\ 2 ) 2 dH 
V ~ 45 du 

These two equations are equivalent to 

ri = 



it 



C x n, 



{i+\z\ 2 y dn 

45 dv 



(46) 



(47) 



with n given by ()22|) . (As was already mentioned, the WCL equation ()32|) is of the 
same form as (|45j) . if H is calculated along the unperturbed orbits.) It is worth noting 
that the equations of motion ()44l) . (|46|) can be formulated in terms of the generalized 
Poisson bracket 



dTdU 
<9q dp 



dT_dU (l + u 2 + v 2 ) 2 
~dp~ 9q + IS 



dTdH 
dv du 



dT dH 

du dv 



(48) 



for any function T of the extended-phase-space coordinates. 
According to (fT§|) . the action along a periodic orbit is 



S po {E) 



T(E) 



1/ • • \ ^r, UV — VU 

- p q - q-p + 25 

2 \+u z +v z 



The functional of the second variations 1Z) )0 ,[r],T\ is still given by with 



dt. 



(49) 



r) = (A 5 i/,p,£) 



_L fx 2^/5 fa 

^ r q ' i + u 2 + v 2 ' 



Sp, 



2VS5u 
1 + u 2 + v 2 



(50) 
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as well as the 2(d + l)-dimensional unit symplectic matrix J and the second variation 
of the classical Hamiltonian (1211) 



1 



. . d 2 n nx d 2 n d 2 n 

A t Aj— h 2\iPj——T h piPj- 



dqidqj 



dqidpj 



dpidpj 



+ §(" 2 + £ 2 ) (-uC 1 -vC 2 + C 3 ) 



+ Vs 



(i + N 2 ) 



dn Bn\ SC 



<9v 



' 9m / 



(51) 



The stability matrices and the Maslov indices are determined by the linearized dynam- 
ics rj = JdH^/dr), as usual 021- 

There is, however, one component of the trace formula which is absent for spinless 
systems. It is the quantum phase correction that emerges in the semiclassical limit 
of the spin-coherent-state path integral [31 j . Originally derived by Solari [HO], Ko- 
chetov [221) an d Vieira and Sacramento [34| for the semiclassical spin propagator, it is 
given by 



SK 



1 f 1 



B(t)dt, 



where 



B{t) 



d {i + \z\ 2 fdn 



Oz 



2S 



dz 



+ c.c. 



(52) 



(53) 



is evaluated along the classical path. Being purely of a kinematic nature, this phase 
plays the role of a normalization [221- For the Hamiltonian 1)211) we find 



B(t) = -uCi - vC 2 + C 3 = 2A(t) - [(f(t))*z + f(t)z*]. 



(54) 



In the semiclassical limit of the partition function (|18|) the Solari-Kochetov phase is 
calculated along the periodic orbits. Since (p SK ~ K°, this phase survives the Fourier- 
Laplace transform from T to E without affecting the stationary-phase condition and 
thus enters the trace formula: 

Sg sc (E) = Ao(E) cos (-\s po (E) + ^ (E) - | ffpo ) . (55) 

po ^ ' 

Note that, as follows from (|54[). (1/2) r](E) in the WCL result ([3"U[) corresponds to 
the semiclassical Solari-Kochetov phase. 



4.4 Applicability conditions and coupling regimes 

We now will show that the semiclassical approach of the previous section is not re- 
stricted to large spin, provided that the Hamiltonian is linear in spin. Thus, we will 
go beyond the formal limit (|43[) and present the arguments justifying this step. We 
also identify the possible regimes of weak, strong, and intermediate coupling that are 
unified within the semiclassical treatment of the extended phase space. 
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The trace formula (|55jl results from a stationary-phase evaluation of the path inte- 
gral (|18j) for the partition function. In the discretized version one finds the stationary- 
phase points in the space of {2d+2)N variables q(fy), p(fy)> / = 1, ■ ■ ■ ,iV [cf. PU|)]. 
Each stationary point corresponds to a periodic orbit that satisfies the semiclassical 
equations of motion (|44j). (|45|1. In this section we evaluate the path integral in a dif- 
ferent, but equivalent, way. Namely, we first integrate over the spin variables z{t{) by 
stationary phase, keeping the orbital part fixed. The result, which consists of an expo- 
nential and a prefactor, depends on q(fy), p(tj)- Then we integrate it over the orbital 
variables, again by stationary phase. At the end we should obtain the same sum over 
the periodic orbits (f53|) , but the two-step calculation allows for a better control of the 
approximations made. 

The described procedure is very similar to what was done in Sec. 14.11 starting 
with the representation (|26|) for Z(T). The difference is that now A4 s [q, p;T] of (|27[) 
is evaluated by stationary phase. A similar calculation was done in Ref. where 
the trace formula for a spin evolving in the external magnetic field C(q(t),p(i)) was 
derived within the approach of the previous section (including the Solari-Kochetov 
phase). The result is valid even if (q(i),p(i)) is a non-classical path, as in the present 
case. Although the final formula in j^Hj was written for spin s = 1/2, it can be 
straightforwardly generalized to arbitrary spin. Transforming it from the energy to the 
time domain, we obtain the stationary-phase result for 7W s [q, p;T] which coincides 
with the quantum-mechanically exact expression (|M|) for any value of spin. 

It is clear now that the stationary-phase evaluation of 

f V[q\V[p] r - r „ 1/fe1 sin{(g + l/2)77[q,p;r]} 
Z(T) = J 72^F ex P {^[q,p;T]/fc} sin{r?[q;p;r]/2} (56) 

yields the trace formula (|55jl . The stationary-phase condition is fulfilled by the classical 
periodic orbits (q(i),p(i)) which satisfy the equations of motion (|44|). We will also 
need the representation for the partition function, in which the orbital part of the path 
integral (|TH|) is evaluated by stationary phase, producing 

Z(T) = j Vfi s [z* , z\ ^ exp {Hip [q po , p po , z*,z;T]/h + isUi [q po , p po , z* , z;T}} 



po 



J?^e W {iK$[r, OTh ,z*,z-T]}, 



(57) 



where the periodic orbits in the (q, p) space are determined by the Hamiltonian (|21j) for 
a given non-classical path z(t). The integral of second variations of the orbital variables 
is defined in (|10 |) -(|12 |) with the full Hamiltonian 1)2 1|) (the vector of small variations 
^7orb ought not to be confused with the functional ??[q, p;T].) Again, applying the 
stationary-phase approximation to (|57|l. we should recover the trace formula (|55|l. 

With the help of ()56|) and (|57|) we can consider the effect of spin-orbit interaction 
on the orbital degrees of freedom. Assuming that for the orbital Hamilton's principal 
function the semiclassical condition TZo/h 3> 1 is always fulfilled and r],TZi ~ |C| T or b, 
we distinguish the following regimes: 

(i) TZo/h S> ST] (WCL). The orbital classical dynamics is determined by TCo, i.e., the 
spin-dependent factor of ()56|) does not influence the stationary-phase condition. Thus, 
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the results of Sec. 14.11 are recovered. If, in addition, |C| Torb 

» 1 along the periodic 

orbit, T) is given by (|39j) and the fast spin precession can be eliminated from the action. 
The same conclusions follow from (|57[). where the periodic orbits will not depend on 
the spin path z(t), and we can pull the periodic orbit sum and all the factors, except 
for exp[is1Zi], out of the path integral. 

(ii) IZo/h ~ rj (SCL). In this case both TZq and 1Z± determine the stationary-phase 
condition in (|57j) . The fluctuation integral does not explicitly contain h and is slowly 
varying. Since the adiabatic condition |C|T or b ~ T^-o/h ^ 1 is satisfied, the orbital 
motion is not sensitive to small variations of the spin path z(t). This means that the 
second- variation integral after the stationary-phase integration over T>[i s [z*,z] will be 
determined only by the spin-dependent part of 1Z\, i.e., 



where all the quantities are taken along the periodic orbits, and the integrals of second 
variations are calculated separately for the orbital and the spin variations. The Solari- 
Kochetov phase r]/2 was added to the exponential. The Fourier-Laplace transform (jHJ) 
of Z SC (T), done by stationary phase, yields the trace formula. The stationary-phase 
condition E = —dlZ/dT seems to be spoiled by the Solari-Kochetov phase. However, 
as shown in Appendix El the periodic orbits form continuous families on the Bloch 
sphere in the adiabatic limit. Therefore, r/, being the angle of rotation of the Bloch 
sphere (cf. Sec. 14.11 and Ref. [BE]); is a multiple of 2ir. Thus, rj remains constant as 
T is varied, i.e., E = —dTZ/dT = —dlZo/dT (this does not mean, of course, that the 
function E(T) is independent of spin, since TZq is evaluated along the periodic orbit 
in the extended phase space). After the trace formula is resummed to eliminate the 
spin part of action and stability denominator (cf. Appendix El, the result (|41|) will be 
reproduced. Clearly, the orbital integral of second variations in (|58f) will be transformed 
to the stability determinants for the periodic orbits of (fl"2)) . 

(iii) 1Zo/fr ~ r], but in some regions of phase space the magnetic field |C(q, p)| 
becomes small, so that the adiabatic condition fails (intermediate coupling). This is 
the regime when the semiclassical theory in the extended phase space is especially 
useful, since both the WCL and SCL formulations are not applicable. In particular, 
the general trace formula l)55j) should remain valid if the orbit contains mode-conversion 
points where C(q, p) = 0. In the weak-coupling regions the influence of spin degrees 
of freedom on the stationary-phase point is negligible, of course. 

(iv) S ~ IZo 3> h (large-spin limit). This is the standard case when the spin-phase- 
space semiclassics generally works, even if the Hamiltonian is nonlinear in spin. The 
orbital motion is affected by spin, and the trace formula (|55|) is appropriate here. 

We conclude that for a Hamiltonian linear in spin, the general semiclassical re- 
quirement of large spin S 3> h, under which the trace formula h55\) was derived, now 
becomes unnecessary. In fact, in many physical problems the spin is small (S ~ H) 
and the Hamiltonian is linear in spin. The trace formula (|55[) can be used in all four 



Z SC (T) = ^2exp{zR, /h + i(s + 1/2)77} 



po 



e W {in^[r] ovh ,z*,z;T]} 




(58) 
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regimes. It can be simplified in the limits of weak and strong coupling, but, since these 
are asymptotic limits, their boundaries may not be strictly defined in specific numer- 
ical examples. The trace formula in its general form ensures that all the effects of 
spin-orbit coupling on the density of states are taken into account in the semiclassical 
approximation . 1 

It can be inferred from the discussion in Appendix El that the number of periodic 
orbits in the extended phase space increases with the increase of spin-orbit interaction 
strength. As a consequence, the number of bifurcations may rise, making our theory, 
in a sense, technically more challenging. However, we believe that the method should 
work reasonably well when the number of relevant orbits is not too high, e.g., in the 
intermediate coupling regime (see the example in Sec. l5.2j) . An extension of our theory 
that includes the treatment of bifurcations by uniform approximations is desirable. 

As was already mentioned at the end of Sec. El there is a freedom in choosing the c- 
number Hamiltonian for the path- integral representation of the partition function (|18|) . 
It is important, therefore, to understand how the semiclassical theory of the previous 
section depends on a particular choice of this Hamiltonian. In our path-integral con- 
struction we assigned to the spin operator hs the classical symbol hsn. The difference 
between this symbol and other possible symbols is 0{h) (not 0(hs)\). Thus, different 
symbols of the spin Hamiltonian H s disagree by 0(fo|C|), which is of the order of 7i s 
itself, unless the spin is large. This means that outside of the large-spin limit the 
results of Sec. 14.31 are valid, generally speaking, only with the specific choice of the 
c- number Hamiltonian (|21[) . Note that even in the large-spin limit, a change of the 
symbol will show up in the phase of the trace formula, since the action will be modified 
by 0(h\C\T). 

5 Applications of semiclassical methods 

We shall now apply the semiclassical methods of Sec. 0] to two specific systems. In 
Sec. 15.11 we study the free two-dimensional electron gas with the Rashba Hamiltonian 
and the Jaynes-Cummings model, that allow analytical treatment on both quantum- 
mechanical and semiclassical levels. Then in Sec. 15.21 we consider a numerical example 
of a quantum dot with harmonic confinement and Rashba interaction. This system 
is a good test case for our new semiclassical approach in the extended phase space, 
since the SCL method suffers from the mode-conversion problem, while the WCL trace 
formula completely neglects the spin-orbit interaction. 

5.1 Rashba and Jaynes-Cummings models 

The free two-dimensional electron gas with a Rashba spin-orbit interaction [8] in a 
homogeneous magnetic field B = i?o e 3 is characterized by the Hamiltonian 

1 At this stage we have to restrict ourselves to the case of C(q, p) depending on cither qj or pj, but not 
both, for each j = 1, . . . , d, due to the problem of the "no-name" term (see Sec. I4.2|l . 
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Using the symmetric gauge for the vector potential, the components of noncanonical 
momentum are given by fri = p\ — (eB§j2c) q 2 and % 2 = p 2 + (eBo/2c) q%. an is 
the Rashba constant [HI^Bj m * * s the effective mass of an electron, e is the absolute 
value of its charge, g* is the effective gyromagnetic ratio, [±b = eh/2mc is the Bohr 
magneton. 

Due to the commutation relation [712, ^l] = iheBo/c we can introduce canonically 
conjugated operators Q = —Tr 2 ^c/eBo and P = —tti ^Jc/cBq satisfying [Q, P] = ih. 
Then, the Hamiltonian (|59|) can be written as 

H = y (Q 2 + P 2 ) + 2Hku c (Qs 1 - Ps 2 ) + %w c s 3 , (60) 

where we have introduced the new parameters 2 k = or \J eBo jlit^fc, 7 = g*m*/2m, 
and the cyclotron frequency uj c = eBo/m*c. 

The Jaynes-Cummings model |19j . that describes the simplest possible interaction 
between a bosonic mode and a two-level system, has a similar Hamiltonian, although all 
its ingredients have a different physical meaning. It can be recovered after subtracting 
huj c / 2 from the Hamiltonian (|60|) and taking s = 1/2. 

For s = 1/2 the quantum- mechanical energy spectrum of (|6U|) is known analyti- 
cally 0: 



E = fno c {l--f)/2, E± = hw c [ k ± V(l - 7)74 + 2khn 2 ) , k = 1,2,3,... (61) 



Using Poisson summation, it can be identically transformed to an exact quantum- 
mechanical trace formula [2]. The smooth part is g(E) = 2/huj c and the oscillating 
part becomes 

, , 2 v-^ / hn 2 \ 



^ c ^ \ ^(1 - 7) 2 /4 + 2Ek 2 /oj c + ft 2 K 4 / 

00 r / e \ 

cos 2irk ( — + /iK 2 ± ^/(l - 7 ) 2 /4 + 2£k 2 /w c + h 2 K 4 J 
fc=i L V c / 



(62) 



One can also find analytically the oscillating part of the level density in the WCL and 
SCL 0. 

To illustrate the semiclassical methods discussed in this paper, we will generalize 
5g^. CL to arbitrary s, first, using the WCL approach of Sec. 14.11 and, second, working 
in the extended phase space (Sec. I4.3j) in the WCL. 

5.1.1 WCL method 

According to (|21|) the phase-space symbol of the Hamiltonian (|6U|) is 

n = y (Q 2 + P 2 ) + 2hsKUJ c (Qni - Pn 2 ) + hsjuj c n 3 . (63) 



2 This definition of k is applicable only within Sec. 15.11 
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It is suitable to introduce the complex variables a = {Q + iP) / V2 and n = ni + iri2 
and to rewrite the Hamiltonian (|63|) in the form 

TC = uj c a*a + V2hsK,uj c (an + a*n*) + hs^uj c n^,. (64) 

In the WCL we can identify two small dimensionless parameters: fko c /E (validity 
of semiclassics) and Hk 2 (weak coupling). We impose that hu c /E ~ Hk 2 <C 1, meaning 
that in the double limit E — > oo, k — > the combination En 2 /oj c is kept constant. It 
corresponds to the formal limit h — » 0. 

We need to solve Eq. (32) with A<® = 7w c /2 and /(°) = k^IEu^ exp(-iu; c i). Its 
periodic solutions are 

z = — |z| exp(— iuj c t), (65) 
where |z| is constant and is given by the quadratic equation 



\z\ 2 - \_L |z|-l = 0. (66) 



Of the two roots 



1 1 2k^/2E^ c V 8Ek 2 y ' 

we will choose |z| + (|z|_ yields the same result) and calculate the integrand of (|31|) 

2A(°) - (/(°>z + / (0) z*) = w c (l + 2A), (68) 

where 

/(I-7) 2 2Ek 2 , . 

V 4 Ld c 

The unperturbed trace formula for the system without spin degrees of freedom, 
corresponding to Hq = 7v 2 /2m* , is that of a one-dimensional harmonic oscillator with 
the cyclotron frequency uj c and reads 6 

^ (0) (^) = ^f;(-ircos(r^). (70) 

Integrating (|68j) over the r repetitions of the primitive period of the unperturbed system 
Tp^ = 27r/u> c , we find 

ri r (E) = 27rr(l + 2A). (71) 
Then, according to (|29|) . the WCL trace formula is 

^ CL (^) = ^- Y(-iy^) Sln[(2s + 1)2 " rA] cos f ^ . (72) 
ysc v 7 ?^cfr{ sin[2vrrA] \ Hu c J 

For spin s = 1/2 it yields 

2 



r=l 

oo 



± r=l 



27rr ( - — ± A 

njJJr 



(73) 



in agreement with the result of 
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5.1.2 WCL in the extended phase space 

Now we will carry out the procedure outlined in Sec. 14.^1 for the Rashba Hamiltonian. 
The implementation of the WCL in the extended phase space for a general Hamiltonian 
of the form (j!7j) can be found in Ref. . 

The classical dynamics of the system is described by the equations of motion (|44j) . 
(|47|) . which now have the form 

a = —iiv c (a + v2HsKn*), 

n = iuj c (^n — 2V2Kn 3 a*), (74) 
ri3 = —iy2KUJ c {na — n*a*). 

Along with the energy conservation TL = E there is another conserved quantity \a\ 2 + 
hsns = const. It means that the system (|63|) is classically integrable. Its general 
analytic solution has been given in pQ. Here, however, we are interested in the periodic 
solutions, which infer the knowledge of certain initial conditions. Let us make an ansatz 
for the particular isolated periodic orbits with n 3 = const: 

a = \a\e- int , n = -£\n\e im (£ = ±1) (75) 

with constant parameters f2, |q| and |n|. Plugging (|75jl into (|74jl we obtain the relations 
between them: 

ft = u c (l-y/2£hsK\n\/\a\), 

Q = w c (7 + 2v / 2^n 3 |a|/|n|) ; (76) 
E = u c (\a\ 2 + hs^ns — 2V2^hsn\a\ \n\). 

For the rth repetition of the orbits (|75|) we can calculate the action 

rT ( E "> faa*-a*a ^ zz* — z*z\ , , , 2 



SUE) = <b dt[ + hs = 2irr(\a\ z + hsn 3 ) + 2vrr/is 

Jo V ^ l \t + \ z \ ) J 

= 2nr + Hs(l - j)n 3 + 2V2£hsK\a\\n\\ +2%rhs (77) 
and the stability angle 



AttklJc H 2 s 2 \n\ 2 2\a\ 2 

p° = r —\ ~i6r + ik ~ 2hsn * + 27rr - (78) 



It seems to be difficult to solve the algebraic equations (fTo]) with respect to f2, \a\, 
and \n\. Moreover, from the numerical search of the periodic orbits in this system we 
know that for rather large k there might appear other periodic orbits with more com- 
plicated shapes. Up to now we have not used the WCL conditions in our calculations. 
In particular, the expressions for the action and stability angle are valid for any values 
of parameters. We gain considerable simplifications in the WCL, that corresponds to 
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the formal expansion in a series of h: \a\ = |a| (0) + h\a\^ + . . ., etc. Then we can find 
that, for instance, 



(0) 



V^M, ^ (0) = w CJ (79) 



1(0) 



£(1 - 7)/2A, |n| w = kV2E/A^U~ c . (80) 

We conclude that, in the leading order, the orbital motion is unaffected by spin and for 
the two periodic orbits we have ri£ =1 (i) = — n^ = _i(t), i.e., the spins are opposite at any 
time. Moreover, we know from the general considerations |3S] that in the WCL only 
two periodic orbits are possible for the Hamiltonian (|64|) . and they are given by ()75|). 
For these orbits we can calculate the action = (2ttE/uj c ) + S t hsrj r (E) [Eq. ([71)1] 

and the stability angle A£ = Ar = rj r (E) using (f77|) and (f751) . respectively. As was 
discussed in Sec. 14.31 the Solari-Kochetov extra phase cp^ , which in the WCL is equal 
to (1/2) £ r] r (E), should be added to the trace formula. Thus, the total phase that 
enters the trace formula is (up to the Maslov indices) 

= ?l + ipf K = + (2s + l)2vrr£ (2A + 1) + 0(h). (81) 

Then we find the oscillating part of the level density 



i 00 t(o) 

« 0L (*> - * E E -C5 - <«> 



where = 2tt/uj c is the unperturbed primitive period. Adding the appropriate 
Maslov indices [SOI c£ = 1 + 2 [£ ry r (£')/27r], where [x] is the largest integer < x, and 
summing up £ = ±1, we can recover formula (|72|). 



5.2 Quantum dot with the Rashba interaction 

We consider a two-dimensional electron gas in a semiconductor heterostructure, later- 
ally confined to a quantum dot by a harmonic potential. We assume that its Hamilto- 
nian 

2m* 2m* 2 2 v ^ v ; 

includes a spin-orbit interaction of Rashba type [S], where 3 k = aji/h 2 . The c-number 
Hamiltonian (|21|) in this case is 

P? m*fi?g? m*Ongn „ , . , 

= 2m* + 2m* + ^ + + 2S <^ ~ ^ (M) 

We will treat this system within our semiclassical approach in the extended phase 
space. As discussed in j2j, the WCL trace formula fails to account for the spin-orbit 
interaction. Indeed, without spin-orbit coupling the only periodic orbits of the system 

3 This definition of n is applicable only within Sec. 15.21 
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are the two self-retracting librations along the principal axes (we assume fii/f^ to be 
irrational). The effective magnetic field C(p) = 2Ke 3 x p changes its sign together 
with p. Hence the spin precession generated by a libration is self-retracting as well, 
making the rotation angle for the Bloch sphere r/(E) vanish. Then the modulation 
factor Mp (E) = 2s + 1 for both orbits is trivial. Thus, the WCL trace formula is 
that of the two-dimensional anisotropic harmonic oscillator [5] multiplied by the spin 
degeneracy factor (2s + 1): 

2s + lv^ {-\) k f J 2irE\ , . 



In the SCL the two pendulating orbits are still present and contain the mode- 
conversion points [2], where the interaction term vanishes. Thus the SCL approach in 
its original form cannot be applied either. Although one can fix it to some extent by 
an ad hoc procedure ^3], it would be only natural for us to resort to our extended- 
phase-space semiclassics, that does not suffer from the mode conversion problem and 
is not restricted to weak or strong spin-orbit coupling. 

The classical dynamics is described by the equations of motion (|44jl and (|46j) that 
now become 

qi = — + T~ l n~ t o = — + 2SKn 2, (86) 

m* 1 + u z + v z m* 
p 2 4Sku P2 

m* 1 + u 2 + v 2 m* 
v = —k(1 + u 2 + v 2 )p2 + 2ku(up2 — vpi) 
px = -m*n 2 iqi , (89) 
p 2 = -m*n 2 q 2 , (90) 
ii = —k(1 + u 2 + v 2 )p\ — 2kv{up2 — vpi). (91) 

Equations l|88j) and (|9Tj) are equivalent to (|47|). or, 

n\ = 2Kpin 3 , h2 = 2Kp 2 n 3 , h 3 = -2k(piUi + p 2 n 2 ). (92) 

For a numerical study, it is convenient to describe the system in dimensionless 
quantities. Let us choose a characteristic frequency Qo and define the dimensionless 
u)\ = Qi/Qo and L02 = ^2/^0 (note that there is a degree of arbitrariness in the choice 
of Qq). Then we can construct the following scaled variables: 



92 = — ~ ; , -2—2- = — - ISnnx, (87) 



m*tt 1 2Sm* . . 

q=q v^' p=p v^m' K=K v^T' (93) 

(E,H) = ^ r (E,H), t = n t, (S,H) = ^(S,n). (94) 
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The Hamiltonian in the scaled coordinates has the form 



jj P1.P2, u\q\ u? 2 q\ n _ vpi - up 2 



With the scaled action ljl9"j) 
S -- 



1 udv — vdu 
-(p • rfq - q • dp + — — 2"—^ 

2 1 + u z + v z 



(95) 



(96) 



the applicability condition of the semiclassical approach now reads S 3> l/2s. 

The equations of motion for the tilded variables are (time derivatives with respect 
to t) 



Qi 
Pi 
hi 



p 1 + kn 2 , q 2 =p 2 -Rn 1 , 
-oj\qx, p 2 = -u\fa, 
2Kpin 3 , h 2 = 2Rp 2 n 3 , h 3 



-2k(pi«i + p 2 n 2 ). 
Note that the system possesses certain discrete symmetries: 



(97) 
(98) 
(99) 



Rl2 


qi - 


* -qi, 


q\ - 


* -92) 


n 3 - 


> -«3, 


t - 


> -i 


(100) 


Ri 


qi - 


* -qi, 


P2 - 


* -P2, 


m - 


-> -m, 


i - 


* -t 


(101) 


R2 


q.2 - 


+ -Q2, 


Pi - 




n 2 - 


■>• -"-2, 


i - 


•> -i 


(102) 



There are trivial solutions to the equations of motion — the pendulating orbits with 
"frozen spin": 

• The pair of orbits pendulating along the q\ axis with spin n 2 = dbl. The 



phase-space coordinates along these orbits are 

qi(t) = ± (y2E + k 2 /wi) sin(^it), 
v(t) = ±l, q 2 (t) = p 2 {t) =u(t) = 0. 



Pi(t) 



=Fk ± y2E + R 2 cos(cJit), 

(103) 



The period of the orbits is Ti = 2it/u)\. The orbits are invariant under R\ 2 and 
-Ri (R2 produces the symmetry partner, i.e., maps + onto — ). 

• Similarly, there are a pair of orbits pendulating in the q 2 direction with spin 
m = ±1. They are given by 



q 2 (t) = ± 2E + R 2 /uj 2 J sm.(u) 2 i), p 2 (t) = ±R± \J 2E + R 2 cos(u 2 t), 
u(t) = ±1, qx(t) = pi(i) = v(t) = 0. (104) 

The period is T 2 = 2tt/uj 2 . The orbits are invariant under R\ 2 and R 2 {R\ 
produces the symmetry partner). 

Other types of periodic orbits can be found from the numerical solution of the 
equations of motion (|97 j) -(|99 j) . In our numerical example we use u\ = 1.56 and u> 2 = 
1.23. For a large range of parameters with < R < 0.75 and E > 8 we find the 
following non-trivial periodic orbits: 
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Figure 1: Periodic orbits (solid line) and D^ 2 (dotted line) for k = 0.67, E = 60, represented 
in different cross-sections of the phase space. Here Q = ^/2S/m*n , P = ^/2Sm*Q, . Note the 
different scales along the axes. In the lower left plot the two orbits are almost indistinguishable. 
For the n(i) time dependence see |27| . 



• Two pairs of orbits and D^ 2 oscillating around in the configuration space, 
with ri2 ~ (Fig. P). The spin is rotating about ri2 axis. The superscripts (±) 
denote the sense of rotation in the subspace (c/i,^)- 

• Two pairs of orbits and D^ 2 oscillating around A^ , with rii ~ (Fig. |2j). 
The spin is rotating about n\ axis. 

For stronger couplings k > 0.75 or smaller energies E < 8, new orbits bifurcate from 
the A and D orbits. Near the bifurcations the trace formula would have to be modified 
by uniform approximations The periods of the orbits are shown in Fig. |21 

The spin-orbit strength k depends on the band structure For example, for an 
InGaAs-InAlAs quantum dot with ~ 100 confined electrons one would obtain a value 
of k ~ 0.25. In order to have the effect of spin on the orbital motion more pronounced, 
we choose k = 0.67 for our numerics. 

The stability determinant in the trace formula can be calculated according to the 
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Figure 2: Same as in Fig. ^for orbits (solid line) and Dy 2 (dotted line). In the upper middle 
plot the two orbits are almost indistinguishable. 



general prescription (Sec. 14.3)) . The second variation of the Hamiltonian l|95|) is 

~2 7.2 , ,2\2 ,,2x2 
-jj(2) Px i Pi 1 ^1 A 1 , ^2 A 2 

n = Y + Y + — + — 

~f ~ ~ 2\ ~ 2\[2uv _ . 

+ K(UP2-VP 1 ){£ +V) + K— 5- 2(^2^ -PKJ 

1 + 1i z + U z 

+ fev ^ 1 I 2 T 2 ~ 1 ~ "V2 2 - 1 friz, (105) 

\l + tr + i> z / \1 + u z + v z J 

where the variables A and p are scaled like q and p, respectively. Numerically solving 
the equations of motion for variations 

2uv ( 2v 2 \ 

Ai = P!-kV2 - =- . -1 )u, (106) 



1 + U 2 + V 2 \ 1 + U 2 + V 2 

(2u 2 \ 2uv 

-— r ^-l)Z + RV2 T — r —v, (107) 
1 + u z + v z J I + u z + v z 

2y/2uv r ( 2u 2 \ , .„ . 

V = - K l + u 2 +v 2 Pl + KV2 [ 1 + u 2 + v 2 - 1 ) P 2 + 2 < U P 2 ~ V PlK> ( 108 ) 
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Figure 3: The periods of the A and orbits as functions of energy for k = 0.67. The periods 
are given in units of l/Oo, the energy unit is 2SQq. For the orbits and A^ the periods are 
independent of energy and equal to T\ ~ 4.02768 and Ti ~ 5.10828, respectively. 



h = 




h = 


2 \ 

— uj 2 M, 


i = 


kV2 (j 



(109) 
(110) 

2v 2 \ 2V2uv . . . . 

1 Pi ~ K l~, op2-2K{Up2-Vpi)U, (111) 



1 + V? + V 2 J 1 + U 2 + V 2 

one determines the reduced monodromy matrix and then finds the stability determinant 
det (M po — I4) of the periodic orbits. 

After calculating the Maslov indices within Sugita's approach |32j and the Solari- 
Kochetov phase by Q52|) . we can compute the oscillating part of the density of states 
using the trace formula (|55j) (Fig. 2J. Note that while the classical dynamics in the 
scaled variables is independent of the value of spin, the density of states will depend 
on S, since the unsealed action 2SS enters the phase of the trace formula. We choose 
the physically meaningful S = h/2 in our example. To ensure the convergence of the 
periodic orbit sum, the density of states was convoluted with a normalized Gaussian, 
exp[— (£'/7) 2 ]/7- v /7r, i.e., it was smoothed out with the energy window ~ 7. With the 
averaging parameter 7 = 0.6, the first repetitions of the 12 primitive periodic orbits 
A and D were sufficient in the trace formula. The semiclassical result for the density 
of states is compared with the quantum-mechanical curve, obtained from a numerical 
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diagonalization of the Hamiltonian (|83[). We observe a rather good agreement between 
the two, especially for the oscillation frequencies. The difference in the amplitudes can 
be explained by the vicinity of the bifurcations in the parameter space. Indeed, the dis- 
parity becomes larger near the avoided bifurcations, where the stabilities are extremely 
small (Fig. EJ). In principle, these energy regions should be treated by a uniform ap- 
proximation [55]. The matter is complicated, however, by the fact that the avoided 
bifurcations are non- generic and of codimension larger than 1. The theory for such 
bifurcations is developed only for two-dimensional systems. Our system is effectively 
three-dimensional. In addition to the elliptic and (inverse) hyperbolic orbits, it also 
has the loxodromic orbits. Thus, the extension of the standard theory of bifurcations 
is difficult and still needs to be developed. 

A more detailed view of 5g(E) is shown in Fig. |SJ For comparison, we added there 
the density of states calculated by the WCL trace formula with spin S = h/2, which, 
in this case, is the trace formula of the two-dimensional anisotropic harmonic oscillator 
without spin-orbit coupling [H] multiplied by the spin degeneracy factor 2 (see above). 
Clearly, it is shifted by phase from the exact density of states. It is worth mentioning 
that the dashed (WCL) curve would come very close to the exact quantum-mechanical 
curve if the former is shifted by —k 2 /2 in the scaled energy E. Note that the action 
of the trivial periodic orbits is Si = (tv/oj{){2E + k 2 ). Thus the shift can be related to 
the perturbative correction to the action in the extended phase space, which is of the 
second order in k. To explain the shift, it would be interesting to develop the second- 
order perturbation theory in the parameter k for the trace formula in the extended 
phase space, similar to that proposed in [7j. (The first-order perturbation theory [S] 
should be equivalent to the WCL trace formula.) Apparently, the shift of the dashed 
curve cannot be justified within the standard WCL approach of With the increase 
of k, the non-integrability of the system will show up, and the shape of the WCL curve 
(describing the integrable system) will start to deviate from the exact density of states. 

6 Conclusion 

We have presented a detailed discussion of a new semiclassical method for systems with 
spin-orbit coupling. The key idea of this approach is to introduce an extended phase 
space of orbital and spin degrees of freedom where the semiclassical dynamics takes 
place. The recipe for the construction of such a phase space, as well as the equations 
of motion there, can be obtained from the path-integral formulation of quantum me- 
chanics using the spin coherent states, with a subsequent stationary-phase evaluation. 

We have classified the possible regimes according to the strength of the spin-orbit 
coupling and the value of the spin. When the Hamiltonian is linear in the spin, our 
method works not only in the standard semiclassical domain of large spin, but also 
for any finite spin. From the path-integral formulation we have directly derived the 
trace formulae in the limiting cases of weak and (under certain restriction) strong spin- 
orbit coupling for arbitrary spin. It was argued that the semiclassical approach in the 
extended phase space recovers the limiting behavior. 

The general method was illustrated with the specific examples relevant, in particu- 
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Figure 4: Upper and middle panels: the stabilities det(M po — I4) of the and D ± orbits. The 
curves come close to the zero level, but do not touch or cross it (avoided bifurcations). Lower panel: 
the oscillating part of density of states for k = 0.67 and S = h/2. QM: the quantum-mechanical 
direct diagonalization. EPS: the semiclassical calculation in the extended phase space, including 
only the first repetitions of the 12 primitive periodic orbits A and D (Gaussian averaging parameter 
7 = 0.6). The energy is measured in units of MIq. 
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Figure 5: A detailed view of the lower panel of Fig. ^Jwith the WCL density of states added. 



lar, to the mesoscopic physics of heterostructures. Our analytical and numerical results 
underscore the importance of the Solari-Kochetov phase correction in the coherent-state 
path integrals and the proper evaluation of the Maslov indices in the extended phase 
space. 

The future work in this direction might include the improvement of the semiclassi- 
cal evaluation of path integrals in the case of strong spin-orbit coupling (restoration of 
the no-name term) ; generalizations of the method treating symmetry breaking and bi- 
furcations by uniform approximations; study of models with the Hamiltonian nonlinear 
in spin and driven by an external system, such as the kicked top ^H] i an d further appli- 
cations to specific systems of interest in spintronics, molecular dynamics, and nuclear 
physics. 
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A Exact quantum-mechanical evaluation of the 
modulation factor (1271) 

A.l Exact calculation of the trace of a spin propagator 

We consider a spin in the external magnetic field C(t) = (2Re f(t), — 21m f(t), 2A(t)) 
described by the Hamiltonian 

H(t) = hs- C(t) = h [2A(t)s 3 + f(t)s+ + (/(*))*$-]. (112) 

Its exact quantum-mechanical propagator between the spin coherent states \z{) and 
(z f \ is HU 



[a*(T) - b*(T) zi + b(T)z* f + a(T) z)zp s 



K s (z}, Zi ;T)- \ 9r \2\8/i , u.|2^ ( 113 ) 



where the coefficients a(T), b(T) are to be found from the equation 

(114) 



d_( a b \ _ .( A f \ fa b 
dt\-b* a* ) ~ 1 \ f* -A ) \ -b* a* 

with the initial conditions a(0) = 1, 6(0) = 0. One can observe that the determinant 
of the matrix 

™ ~ ( 4% $% ) <-> 

remains constant: |a(i)| 2 + |6(t)| 2 = 1- Hence, U(t) belongs to the group SU(2) and 
can be represented in the form 

?7(t) = cos^+i(m(t)-cr)sin^, (116) 

where m 2 = 1 and a is a vector of Pauli matrices. Comparing ()115l) and ()116|) we 
deduce that 

n(t) n(t) n(t) 

a(t) = cos-^- +im 3 (t) sin-y^, b(t) = i[mi(t) - im 2 (t)] sin -y-. (117) 



Since 



we can rewrite (|114|) as 



1 r( t) rr-f A ® *® \ (MSI 
-C(t)-a-[ ), (118) 



±U(t) = - l -(C(t)-a)U(t), U(0)=I 2 . (119) 

Note that this equation has the same form as the equation for a spin-1/2 propagator 
in the standard |s,m,) basis. 
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The trace of the spin propagator (|113|) is 

Z S (T) = TvK s = J dns{z'\z')K s {z} = z'*, Zi = z';T), 

where dfi s (z'* , z') is given by (|16() and 

~a*(T) - b*(T)z f + b(T)z'* + a(T)\z'\ 2 ^ 2 " 



K s {z'\z'-T) 

After the stereographic projection [cf. ((221) ] 

2z'* 



1 + U'l 2 



tl x + in' 2 = sin 9' exp(i(j)') 



n' 3 = cos 9' 



1 - \z 



/|2 



l + |z'| 2 ' 6 l + |z'| 2 

and simple transformations we can cast (|121() into the form 



K s (n';T) = cos 



r,(T) 



Since the measure 



^ , in ■ m(i ) sin — - — I 



2s + 1 

dfi s (z'* , z') = dfj, s (n) = — s'm 9' dO' d(p' 



2 s 



is invariant under rotations, we can choose the axis e3 along m(T). Then 



Making the following transformations in ([12 

2s + 1 



Zs{T) 



2s 



2s + 1 



7i(T) , r?(T)Y 

cos — h z cos 0' sin — - — j sin 9 dO' 

ry(T) . V (T)\ 2s , 

cos h ia; sm dx 

2 2 J 



1 



r?(T) . . r/(T)\ 

- . cos h «x sm 

2isin[r?(T)/2] V 2 2 / 



2s+l 



x=-l 



we obtain the final result 



Zs{T) 



sin[(s + l/2)7y(T)] 
sin[77(T)/2] 



(120) 



:i21) 



(122) 



(123) 



(124) 



Z S (T) = J dfi s (n')K s (n';T) = J dfi s (n') (cos ^ + in' 3 sin 2^ . (125) 



(126) 



where r/(T) is found from 2 cos[t?(T)/2] = 2Re[a(T)] = Tr[f7(T)]. 

In order to give the geometric interpretation of r)(T), consider an equation for the 
spin s propagator [cf. ((119)) ] 



jV s = -i(C-3 s )V s , V s (0) = I 2s+1 , 



(127) 
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where J s = (</f ; </| > "^D are the generators of the (2s + l)-dimensional irreducible 
representation of SU(2). Since V S (T) = exp[im(T) • J s r)(T)], as well as U(T) = 
exp[fm(T) • crr](T)/2], rj(T) has the meaning of the rotation angle around the m(T) 
axis. 

Note that one can alternatively find Z S (T) as Tr[D s (T)]. Choosing m(T) as a quan- 
tization axis and making simple calculations, one would obtain the result coinciding 
with (TTSHl) . 



A.2 Calculation of 77 (T) 

For the reader's convenience we derive ()31|) within our notation. We will closely follow 
the discussion in 5,, where this result previously appeared. 

Without loss of generality we can choose the basis in which e3 = m(T). One can 
decompose U(t) (|115JI at every t into the matrix product 



U(t) 



cos(0/2)e-^ -sin(0/2) 
sin(0/2) cos(0/2)e^ 







(128) 



where a = cos(#/2) e l ^~^ and b = — sin(#/2) e~ % ^ . Decomposition (jl28|) corresponds 
to the choice of a certain section in the principal U(l) Hopf bundle over § 2 , ip being a 
fiber coordinate. Imposing the time periodicity on Ui(t) and recalling that U(0) = I2, 
we establish that 

^(T) - V(0) = r/(T)/2. (129) 



Note that the decomposition (|128j) is not well defined at t = 0, i.e., when 9 = 0. Hence, 
the initial value ^(0) is not determined. (Still, rj(T) can be found unambiguously.) 
Although not essential for our purposes, the problem can be fixed if we choose a different 
decomposition near 9 = 0, related to the former by a gauge transformation [5j I35| . 
We define the projection of the Hopf bundle (9, <j>, ip) S S 3 — ► n € S 2 by 



n • a 



Ua 3 U^ = U^UI 



which is equivalent to n\ + in-i = sin#exp(i(/>), 723 = cos I 
n(t) and ijj(t) from (|119f) . Thus, one should calculate 



(130) 

One can find equations for 



n 



2 



to recover the equation 

n = Cxn, n(0)=n(T). 
After the stereographic projection (|2*2*|) it becomes 

z = -i[2Az + f-f*i?], z{0) = z{T). 

To derive an equation for ip one should calculate 



i> - -(1 + cos 



-*Tr 

2 



(131) 
(132) 
(133) 

(134) 
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Figure 6: Classical spin precessing about the instantaneous magnetic field. 



taking into account (|128|) and (j!33j) . In the language of the fiber bundle theory, it 
corresponds to the pull-back of the canonical connection form [r.h.s. of ()134|) ] onto the 
section defined by (|128|) . The result of the calculation is 

1 77 — 7 7 1 

^ = 27 ITnF " 2 C ' n(z *' z) = [A ~ {f * z + fz * )/2] - (135) 

Integrating it we find 

V (T) = 2[V(T)-^(0)]= f dt 

Jo 

= [ dt[2A-(f*z + fz*)}, (136) 
Jo 

where z(t) and z*(t) are to be found from . 

B Adiabatic limit in the extended phase space 

Here we present an intuitive derivation of the trace formula (|41j) within the framework of 
our semiclassical theory in the extended phase space fSec. l4.3|) . In the adiabatic regime 
the classical system, described by the Hamiltonian H (|21j) . possesses an adiabatic 
invariant — the angle [5 between the instantaneous magnetic field C and the classical 
spin vector n. To demonstrate this, let us consider the quantity a = cos/? = n • nc, 
where nc = C/|C| (Fig. |EJ). Its time derivative averaged over the fast spin motion is 

(d) = (n ■ nc) » (137) 

taking into account that ri, he -L nc- 

Having an approximate integral of motion a has two consequences. First, the fast 
spin precession about C can be separated from the orbital motion in the adiabatic 
limit. The latter then is described by the polarized Hamiltonian 

W Q (q, p) = Wo(q, P) + MC(q, p)| (138) 



zz 



z~z 



ill + z 



C-n(z*,z) 
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with a as a parameter. Second, in a system with one integral of motion in addition 
to the energy, the periodic orbits are not isolated, but appear in one-parameter fam- 
ilies instead. Hence, we should use the modified trace formula jlUj that involves the 
summation over the periodic orbit families (p.o.f.): 



5g sc (E) = -Im^- 



iirh 

p.o.f. 



2vr 



(2mh) 1 /' 2 





1/2 




dJ 







T ppo exp (|«S + iip bK - i\o 
73 x 1/2 

det [M orh - I 2 (d-i) 



(139) 



where the factor in the square brackets is absent in the standard Gutzwiller trace 
formula for isolated orbits, the stability matrix M or b is calculated only for the orbital 
subspace, the action S, Solari-Kochetov phase ip SK , and Maslov index a are evaluated 
for any representative of a family. J = hsa is the classical-spin projection on the 
direction of the field and 

T7(a) = »7[q,p;r£>] (140) 

is the total angle of rotation of the Bloch sphere for a periodic orbit of TC a , where the 
r.h.s. of 1)140(1 is calculated according to (|36]> (cf. Sec. 14.11 and Ref. 36 ). Obviously, 
we have rj 3> 2tt, since the spin vector makes many rotations during the period of the 
orbital motion. 

As a result of the approximate separation of variables, to any classical trajectory 
of H Q (q, p) = E corresponds a family of classical trajectories in the extended phase 
space. This family will consist of periodic orbits if the respective orbit of 7i a is periodic 
and the Bloch sphere makes an integer number of rotations during the orbital period. 
The latter condition follows from the fact that the whole family, not just an isolated 
orbit, is periodic on the Bloch sphere. Thus, we have a requirement 

n(a) = 2vriV(a) (N(a) > 1) (141) 

with an integer N(a). Clearly, this condition can be satisfied only for selected values of 
a. Now we will show that the allowed values of a form discrete sets that become very 
dense in the adiabatic limit. Consider a periodic orbit of T-L a =-x- Let us assume that 
as a changes from —1 to 1 for fixed energy, the orbit can be continuously deformed 
while maintaining its periodicity and that rj(a) increases with a. If for a = — 1 the 
requirement ()141() is not fulfilled, we shift a by a small amount to a = oq, such that 
7/(ao) = 2ttN. The greater the frequency of precession |C|, the smaller the shift. 
Changing a further to ol\ > cxq, we achieve that the Bloch sphere makes one more 
rotation for the newly deformed orbit, i.e., 77(0:1) = 2ir(N + 1). We continue this 
process till we reach the upper bound oil < 1 for which r}(a£) = 2ir(N + L). Thus, we 
have constructed a set {—1 < ctQ < ot\ < ... < ol < 1} of the values of a for which 
the families of orbits in the extended phase space are periodic. We will call it the 
periodic-orbit set (p.o.s.). The better the adiabaticity condition (|3*3*|) is fulfilled, the 
larger the number of elements of the p.o.s. L ^> 1. The periodic orbits within the set 
are related by an almost continuous transformation in the (q, p) space (Fig. EJ). If we 
start with a different periodic orbit of 7i a= -i, we will obtain another p.o.s. The sum 
over the p.o.f. in (|139|) can be represented as the sum over all sets and the elements 
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Figure 7: Elements of the periodic orbit set in the orbital subspace. 



within each set, i.e., 

N+L oo r N+L 

£(•••)=££(•••)=£ £ / dne^(...). (142) 

p.o.f. p.o.s. n=N p.o.s. M=— oo 

In the second equality we applied the Poisson summation formula to the inner sum. 

Our conclusions may seem to disagree with the results of Sec. 14.21 where we have 
found for the periodic orbits in the extended phase space that a ~ ±1. There is no 
contradiction, however, if we recall that in (|37|) the magnetic field C was predetermined, 
so that, in general, the orbital frequency and the overall frequency of precession were 
incommensurate. Therefore the spin must have been aligned with the field in order 
to have a periodicity in the extended phase space. In the present case, on the other 
hand, there is a feedback from the spin to the orbital degrees of freedom, and for 
certain initial conditions the frequencies become commensurate, thus producing the 
periodic-orbit sets. 

We now turn to the calculation of the action which can be divided into the orbital 
and spin parts [cf. (|49j)] . The orbital part S^, is evaluated in the usual way with the 
Hamiltonian TC a (q, p) = E. The spin part is given in terms of the spin orientation 
n(M)by 

S ( ^} n = hs j P ° (1 + cos0) 4>dt = hs[rj(a) + a$(a)], (143) 

where $(a) = j> po |C(q, p)\dt. There is a direct analogy between the problem of 
spin precession and the problem of a spinning top with fixed base point whose axis 
of rotation is moved along a closed circuit. Namely, the magnetic field orientation 
n c*(#c> 0c) corresponds to the axis of rotation, its magnitude |C| becomes the angular 
velocity of rotation, and the spin orientation n(6,(p) is a direction on the top making 
angle (3 with the axis (Fig. [BJ. Then the standard theory of adiabatic motion JJj 
provides the relation 

v ( a ) - 4vr = $(a) - Q(a), (144) 

rp(a) 

where Q(a) = j> po (1 + cos 9c) (pcdt is the solid angle swept over by the axis of 
rotation in the lower hemisphere during the period. Using (|144j) in (|143|) and taking 
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a = cos (3 




Figure 8: Illustration of the coordinates' notation. = J Q * |C(q, p)\dt. 



into account (|141|) . we obtain the spin action 



:(«) 



inks 



. S|Ull 2vrfe(a + l)iV(a) + feaO(a), (145) 

where N(a) was shifted by —2. The second term on the right is the origin of the Berry 
phase. 

The Solari-Kochetov phase ip^ K = r/(a)/2 = irN(a) is equal to an integer multiple 
of 7T. We expect that together with the Maslov index a, it will yield, at the end of 
the calculation, the usual Maslov index for the trace formula of the Hamiltonian 
Ho,(q, p) in the orbital subspace. However, at this point we have not succeeded in 
establishing this connection. 

We can do the summation over the p.o.f. in the trace formula (|139[) with the help 
of (|142jl . where the integral can be computed by stationary phase. For this we express 
all a-dependent quantities in terms of N(a) and then substitute it by a continuous 
variable n. Only sf^l and 2irhs(a + l)N(a) determine the stationary-phase point, the 
remaining functions are slowly varying with n. The stationary-phase condition is 

= 2ttM + -— [S^l + 2irhs[a{n) + l]n) = 2vrM + 2yrs[a(n) + 1], 
where the result of the classical perturbation theory 



(146) 



OS, 



orb 



-dt = -hs$(a) = -2nhsN{a) + 0{N U 



(147) 



da J da 

was employed. Thus, the equation for the stationary point a(n s t) = —M/s — 1 has a 
solution only if the the r.h.s. does not exceed 1 by magnitude. This makes the M sum 
finite, reducing it to the sum over quantized polarizations: 

-N+L 



At 



I ' dne^ iMn+i ^ +2nhs[a{n)+1]n ) {...) 



E 



g ft "orb 



{2mh) 1 ' 2 



dJ 



1/2" 



2vr 



(...) 



(148) 



a=m s /s 
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m s = —M — s is either integer or half- integer, depending on the spin s. It labels the 
quantized spin projection on the np axis. The prefactor due to the stationary-phase 
integration, that appears in the square brackets in (|148|) . cancels the respective factor 
in Q139I) . The outer sum over the p.o.s. in (|142[) . if exchanged with the sum over m s , 

becomes the sum over the periodic orbits of polarized Hamiltonians W Q=ms / s (q, p) = 
,{m s )( 

formulae: 



Wiff s) (q,p) Thus the density of states is given by the sum of the polarized trace 



1 8 

— y y 



rp(m s ) 

-L ppo COS 



-s po 



2 a P° 



(m a ) 



det M 



r(m s ) 



1 po 



l 2(d-l) 



1/2 



(149) 



with the Berry phase 



iff a) = saQ(a) 



= m s 

OL=m s /s jq 



rp{m a ) 
± po 



[l + COS I 



•c 



<ficdt, 



(150) 



which is equivalent to (|4U|) , 

Far from being mathematically rigorous, this derivation illustrates how the sum 
over quantized polarizations can be obtained from the classical-spin dynamics, where 
any polarization is allowed. As in Sec. 14.21 the problem of the missing no-name term 
|25| remains here. A careful evaluation of the Maslov indices is also required. 
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